Scaling Behaviors and Novel Creep Motion of Flux Lines under AC Driving 
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We performed Langevin dynamics simulations for the ac driven flux lines in a type II supercon- 
ductor with random point-like pinning centers. Scaling properties of flux-line velocity with respect 
to instantaneous driving force of small frequency and around the critical dc depinning force are 
revealed successfully, which provides precise estimates on dynamic critical exponents. From the 
scaling function we derive a creep law associated with the activation by the regular shaking. The 
effective energy barrier vanishes at the critical dc depinning point in a square-root way when the 
instantaneous driving force increases. The frequency plays a similar role of temperature in con- 
ventional creep motions, but in a nontrivial way governed by the critical exponents. We have also 
performed systematic finite-size scaling analysis for flux-line velocity in transient processes with dc 
driving, which provide estimates on critical exponents in good agreement with those derived with 
ac driving. The scaling law is checked successfully. 



PACS numbers: 74.50.-|-r, 74.25.Gz, 85.25.Cp 

Introduction. - Dynaniics of elastic manifolds driven 
through a random medium is one of the rich paradigms 
in condensed matter physics. Important examples range 
from ferromagnetic or ferroelectric domains, charge- 
density waves, frontiers of fluids in porous media, Wiener 
crystal, and flux lines in type II superconductors [l|-l3|. 
The notion that the depinning transition by a dc driving 
at zero temperature in these systems exhibits properties 
similar to the critical phenomena of a second-order phase 
transition at thermal equilibrium has been verified suc- 
cessfully in many systems in this class while a pos- 
sibility was raised recently on the absence of divergent 
correlation length below the depinning point [H] . 

Special attentions have been paid to the dynamic phe- 
nomena of flux lines in bulk superconductors with ran- 
domly distributed defects @. It is because, first of all, 
the system is important for applications of supercon- 
ductivity, such as transport of large current and gen- 
eration of strong magnetic field. Flux lines are also 
unique with respect to the depinning transition due to 
the two-dimensional (2D) dynamic degrees of freedom 
under a driving transverse to the direction of magnetic 
field, which hampers an analytic treatment based on the 
functional renormalization group [Tf . Computer simula- 
tions based on Langevin dynamics |8l-[l4| were performed 
in order to fill the gap. 

While the critical depinning properties of fiux lines un- 
der dc driving were addressed quite satisfactorily (see for 
example [T^). behaviors under ac driving have not yet 
been fully explored so far, despite of its importance and 
convenience as an experimental probe. Fortunately, a 
scaling theory for small frequencies around the critical 
dc depinning point has been formulated based on analy- 
sis on an elastic string embedded in 2D random medium 
[l5| , which provides a facet into this problem and motives 
our present study. 

We performed computer simulations for the ac driven 
flux lines in presence of quenched random point-like pins 



based on the Langevin dynamics. Scaling properties of 
flux-line velocity with respect to instantaneous driving 
force of small frequency and around the critical dc depin- 
ning force are revealed successfully based on the scaling 
theory (T5| , which provides precise estimates on dynamic 
critical exponents. From the scaling function we found a 
creep law associated with the activation by the regular 
shaking, with the energy barrier vanishing at the crit- 
ical dc depinning point in a square-root way when the 
instantaneous driving force increases. The frequency ap- 
pears in the creep law, instead of temperature in a simi- 
lar creep motion under dc driving Il6j , in a nontrivial 
way governed by the critical exponents. In order to check 
possible effects from finite systems in simulations on ac- 
curacy of critical exponents, we also simulated two tran- 
sient dynamic processes under dc driving and performed 
systematic finite-size scaling analysis. It is found that 
the estimates on the critical exponents converge nicely, 
and moreover satisfy the scaling law [l7j . 

Model and simulation technique. - The model system is 
a stack of superconducting layers with a j^erpendicular 
magnetic field, same as that used in Ref. [14|. The over- 
damped equation of motion of the ith vortex at position 
Ti is 

= - ^ V,;7^V(j.^^.)_^ v,(7^P(r,p)+F + Fth, (1) 

P 

with rj the viscosity coefficient. The intraplane vortex 
repulsion is given by the modified Bessel function, the 
interplane vortex attraction is approximated by a spring 
potential between two vortices belonging to the same fiux 
line and sitting on adjacent planes, and pinning centers 
are modeled by Gaussian potentials with typical size Rp 
distributed at random positions but with constant pin- 
ning range and dimensionless pinning strength a; F is 
the Lorentz force uniform over the system but varies 
with time as Fit) = Asmtut, and Fth is the thermal 
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noise force. In this work, the units for length, energy, 
temperature, force, time, and vefocity are taken as Xab, 
deo, deo/kB, deo/Xab, ??A^j/(ieo(= tq), and deo/r]Xab, re- 
spectively, where eo = ^oZ-^^/^o-^ab: ^ab the magnetic 
penetration depth, and d the thickness of the supercon- 
ducting layer. The results will be given in dimension- 
less units hereafter. The dynamic equation ^ is inte- 
grated with the 2nd order Rugen-Kutta algorithm with 
At = 0.01 ~ 0.02. More details about our model system 
can be found in Ref. [T^ . 

In simulations there are totally = 180 flux lines in 
the system of lateral size L x i = 30 x 30 and 20 layers 
in the c axis, with periodic boundary conditions. Intra- 
plane vortex-vortex repulsions are cut off at Tcut = 6. 
Each layer contains A^p = 900 randomly distributed point 
pins with Rp = 0.22 and a = 0.2, associated with a 



vortex glass at thermal equilibrium 2M- -P'or different 
system sizes, the densities of the flux lines and the point 
pins are fixed. In the present work we concentrate on 
zero temperature. 

ac driving. - We first investigate motions of flux lines 
under an ac driving. Double hysteresis loops [l5| are 
observed as displayed in the inset of Fig. 1(a) for u) = 
O.OItt. It is clear that for a finite u the velocity of the 
system is nonzero even when the instantaneous driving 
force is below the critical dc dcpinning force Fcq derived 
from motions of flux lines under dc driving |14| . We 
therefore refer to Fco simply as critical depinning force 
hereafter. 

As the frequency uj approaches zero, the double loops 
shrink to the depinning curve of the dc driving |1J], as 
seen in the main panel of Fig. 1(a) where the instanta- 
neous force dependence of the velocity is shown around 
Fco in the increasing branch. According to the scaling 
theory [l^, the behavior of the system for w < cjp = 
Fcq/v^ — O.OStt is governed by the critical properties of 
the depinning fixed point 
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FIG. 1: (color online), (a) Hysteretic loops of the velocity of 
flux lines under ac driving F{t) = 0.4 sin(0.0l7r/:) (inset) and 
the part around Fco in the increasing branch (dashed rectan- 
gular in the inset) for several frequencies (main panel), (b) 
Scaling plot according to Eq.(l2]) (main panel) and asymptote 
of the scaling function (i>{x) in Eq.Q as a; — s- — oo (inset). 
Simulations were performed with L = 30 and numbers of 
samples are 80-200 depending on the frequency. 



here the critical exponents are defined by the onset of the 
velocity v ~ {F/ F^q — 1)'^ under dc driving, the divergent 
correlation length ^ ~ |F/i^cO — and the growth of 

correlation length with time ^ ~ i^/'^ at F ^ Fco. It is 
requested that (f>{x) ~ for x — > cxd in order to recover 
the steady depinning behavior. 

The scaling plot according to Eci.® sees a good col- 
lapsing with l3/i^z = 0.40 ± 0.01 and l/i^z = 0.53 ± 0.02, 
and Fco = 0.232 ±0.002 as displayed in the main panel of 
Fig. 1(b). The estimates on the exponent f3 = 0.76±0.02 
and Fco agree well with the previous results based on dc 
driving at low but finite temperatures {(3 = 0.75 ± 0.01 
and F^o = 0.232 ± 0.001) jH. 

A clear asymptotic behavior is observed for the scaling 
function (j){x) as x — > — oo in the main panel of Fig. 1(b). 
As replotted in the inset of Fig. 1(b), the asymptote is 
well described by 



v{t) 



uj'^/"' exp 



0.35v/l - F{t)/F,o 



.l/2i/z 



(3) 



for F{t) < Fco- The motion of fiux lines is creep like, with 
an effective energy barrier 0.35 y/l~F{t)/Fco which van- 
ishes when the instantaneous driving force approaches 
Fco from below. The creep motions of fiux lines here are 
caused by the regular shaking of the ac driving, instead 
of the thermal activations at finite temperature [3, ■ 
Because of the random pinning potential landscape, a 
regular shaking contributes to activation of flux lines over 
pinning barriers in a complex way governed by the criti- 
cal properties, both steady and dynamic, of the system, 
which is captured by the frequency dependence in the 
creep law A theory on the square-root suppression 
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FIG. 2: (color online), (a) Traveling distance of the flux lines 
after relaxing from a perfect triangular lattice for several flnite 
systems (main panel) and scaling plot based on the flnite- 
size scaling ansatz ((S)) (inset), (b) Asymptote of the scaling 
function S{x) in the scaling theory ((5]) for x — >■ — oo and (c) 
for X Xc- Numbers of samples are 10-50 depending on the 
system size. Error bars are given for the system of L = 30 
and those for other system sizes are similar. 



of the energy barrier is not available at this moment. 

To end this part we notice that, as seen in the scaling 
theory ([2]), the critical exponents v and z appear in the 
form of product in the dynamics under ac driving, which 
leaves a full description of the critical dynamics unavail- 
able. On the other hand, because of the hysteretic re- 
sponse of flux lines to ac driving, it is still hard to make 
any serious analysis on finite-size effects on the accuracy 
of the critical exponents with the computing resources 
available at this moment. These issues are addressed in 
what follows. 

Onset of collective pinning. - Let us now investigate a 
process associated with the onset of collective pinning, 
which involves only the exponent v. We lay at i = 
a perfect triangular lattice of flux lines on the random 
medium, and start to drive the flux lines by a do force. 



The flux-line lattice deforms during traveling when the 
individual flux lines adapt to the random potential land- 
scape. As the result, the dragging caused by random 
pinnings becomes stronger due to the collective pinning 
mechanism jlSj. When the dc driving force is below the 
critical depinning one, the flux lines should be stopped as 
a whole after traveling a certain distance. This traveling 
distance depends on how far the driving force is from the 
critical point, and is determined by the critical properties 
of the depinning fixed point. 

We have simulated this process in several finite sys- 
tems, and the results are displayed in the main panel 
of Fig. 2(a). It is clear that, for a given finite system, 
the traveling distance diverges (since a periodic bound- 
ary condition is adopted) as the driving force approaches 
a critical value given by the system size. The critical 
force decreases as the system size increases. The critical 
behavior for an infinite system is given by 



D^{1- F/F^oY 



(4) 



with an exponent and the convergence from finite 
systems to the infinite system should be described by the 
finite-size scaling theory fl^: 



D = L^^^^S \iF/Fco - 1)1^^" 



(5) 



The scaling plot according to Eq. ([5]) is shown in the inset 
of Fig. 2(a). From the successful scaling plot, we obtain 
lyo/i' = 1.00 ±0.05, f = 1.06 ±0.04, and F^o = 0.2365 ± 
0.0005. The present estimate on the critical depinning 
force agrees with the above one derived for ac driving 
and that in the previous work for dc driving fl^l , but with 
the highest precision since finite-size effects are taken into 
account here. 

Two features of the scaling function S{x) are observed 
as follows: (i) As shown in Fig. 2(b), ^(a;) ~ (— a;)~''° 
for a; — >■ — oo as requested by the scaling theory, which 
achieves the critical behavior (ii) As shown in 

Fig. 2(c), S{x) ~ a{xc — x)^^ for x — > Xc, with a ~ 0.53, 
Xc — 0.89, and p ~ 1.86. This relation gives the system- 
size dependence of the critical depinning force, and the 
divergence of the traveling distance in any flnite system 
with an exponent p, different from that in an inflnite 
system. 

The successful finite-size scaling plots in Fig. 2 includ- 
ing the exponent imply the divergent correlation length 
both above and below the critical depinning force with 
the same exponent. In contrast, absence of diverging 
correlation length below the critical depinning force was 
discussed for an elastic string embedded in a 2D random 
potential landscape The reason for this discrepancy 
remains as an interesting problem for future work. 

Critical slow down. - As a second transient process as- 
sociated with dc driving, we study the critical slowing 
down of flux lines, which provides a way to estimate 
separately the exponent z. Here a steady state of flux 
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lines under driving force F = 0.6 ^ Fco is prepared. At 
t = the driving force is reduced to F^o- The vefocity 
of the flux lines then decreases with time according to 
V ~ t-P/^^ [2^. We have simulated the critical slow- 
ing down for several finite systems, and the results are 
depicted in Fig. 3(a). The velocity decreases quickly in 
a small system since the corresponding critical force is 
large, as seen in the onset process of collective pinning in 
Fig. 2(a). The finite-size scaling behavior at F = Fco is 
described by [13, [13 



v{t,L) = b^/^vibH^bL). 



(6) 
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The scaling plot is performed successfully as shown in 
Fig. 3(b), and we obtain /3/iy = 0.72 ± 0.02 and z = 
1.80 ± 0.03. The slope of the scaling function at small 
argument is given by P/uz ~ 0.40 as required by the 
scaling theory on short-time dynamics [13, [HI . 

Discussions.- The critical exponents estimated in the 
present study satisfy the scaling law /3 -I- i^(2 — z) — 1 
[iTj . The roughness exponent C can be evaluated by the 
scaling relation i/ = 1/(2 - C) U3 as C = 1-06 ± 0.04. 

In conclusion, with the aid of computer simulation 
and scaling theory, a comprehensive picture has been ob- 
tained for the depinning transition of current driven fiux 
lines in type II superconductors with random point-like 
pins. 



FIG. 3: (color online), (a) Critical slowing down simulated in 
several finite systems, (b) Finite-size scaling plot according 
to Eq.©. Numbers of samples are 900-1200 depending the 
system size. 
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